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Abstract. Lemaitre-Tolman-Bondi inhomogeneous spacetimes can be used as a 
cosmological model to account for the type la supernova data. However, such models 
also give rise to large velocities of galaxy clusters with respect to the cosmic microwave 
background. These velocities can be measured using the kinematic Sunyaev-Zeldovich 
effect. This paper presents a calculation of galaxy cluster velocities as a function of 
redshift for such a model. 
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1. Introduction 

Observations of type la sup er novae [U, [2], along witli tlie usual assumption that the 
universe is homogeneous, indicate that the expansion of the universe is accelerating. 
This acceleration necessitates the presence of an exotic form of matter called dark energy, 
which can most easily be accounted for by a cosmological constant. However, because 
the scale of the cosmological constant is so small (~ lO^^^*' in Planck units) it is natural 
to look for alternative explanations for the data. One such alternative explanation is to 
account for the supernova data using a comsmological model that is inhomogeneous.[3|,|l] 
The simplest such models are the Lemaitre-Tolman-Bondi(LTB) spacetimes. [5l El [7] 
These are spherically symmetric spacetimes with pressureless fluid (dust) as the matter. 
Models of this form can be constructed that account for the supernova data as well 
as the standard homogeneous "concordance" cosmology does. [HI El [THl [HI [12] Thus to 
test the difference between concordance cosmology and inhomogeneous cosmology, a 
different test is needed. Such a test is provided by the kinematic Sunyaev-Zeldovich 
(kSZ) effect. [13j This effect comes about when cosmic microwave background (CMB) 
photons scatter off the gas in a galaxy cluster. The spectrum of the scattered photons 
depends on the velocity of the galaxy cluster with respect to the microwave background, 
and thus an analysis of the CMB can be used to find this velocity. Thus, each LTB 
model must be compared not only to the supernova data but also to the kSZ data, and 
only those models that match both sets of data can be considered viable. 

Analyses of the kSZ effect for particular LTB models have been performed [T^ [T5] 
as well as a general analysis for slightly inhomogeneous models. [16] Unfortunately, each 
LTB model is different, and so each must be seperately subjected to the kSZ test (as 
well as the supernova test) to see if that model is viable. In this paper we will present 
a general method for finding the velocity of galaxy clusters in LTB models. We will 
then apply this method to the models treated in [12]. Section 2 presents the general 
method of calculating these velocities. Section 3 applies this method to the models 
of [12]. Conclusions are presented in Section 4. 

2. calculating velocities 

The metric in the LTB spacetimes takes the form 

ds^ = -dt" + Y^jdf^ + r\de^ + sin2^rf02^ (1) 

Here the area radius r is not one of the coordinates but is instead a function of the 
coordinates t and f. The function / is a function of f. An overdot denotes derivative 
with respect to t while a prime denotes derivative with respect to f. The LTB spacetime 
has a prefered timelike vector u"" = (d/dt)"' which is the four- velocity of the dust that 
makes up the matter in the spacetime. However, an LTB spacetime with a cosmic 
microwave background has another prefered vector, which we will call t", the four- 
velocity of the CMB. Since the spacetime is spherically symmetric, there is also a 
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preferred unit spacelike vector r" such that + is an outgoing null vector and t°- — r°- 
is an ingoing null vector. Expressed in terms of and r", the LTB four-velocity u"" 
takes the form 

M« = 7(t" + vr"") (2) 

where v is the velocity of the LTB observer with respect to the CMB and where 
7 = 1/ Vl — f ^. Now let n"_ and n*]^ be respectively ingoing and outgoing CMB photons 
each with the same frequency uo as measured in the CMB rest frame. Then it follows 
that these vectors take the form 

= cj(t'^ - r") (3) 
n\=uo{e + r'') (4) 

Now define the quantity a by 

a = — (5) 

Then it follows that 

V = -Y—- (6) 

We restrict attention to times when the CMB has decoupled from the matter, so that 
are null geodesies. Thus, to find the velocity of the galaxies with respect to the 
CMB, we need to find the behavior of null geodesies in LTB spacetimes. 

To find the behavior of null geodesies, we first recall some facts about the LTB 
spacetimes. From the Einstein field equation, it follows that 

r' = f + - (7) 

r 

where the function F depends only on r. The density is given by 
F' 

P = (8) 
It is helpful to introduce the quantities a, A and B by 

r = ar (9) 
/ = Af' (10) 
F = Br^ (11) 

Then equation ([7]) becomes 

c? = A+- (12) 
a 

whose solution is 

^-^0= / (13) 

Here to is a function of f whose meaning is the time at which the shell of dust with 
coordinate r shrinks to zero radius. To have a genuine big bang singularity, we chose 
to = 0. We use the coordinate freedom to change f to any function of f to set B to the 
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constant 4/9. With this choice, the particular LTB model is then determined by the 
choice of the function A{f). In particular, when A = the LTB solution becomes the 
standard dust FRW cosmology with a = t^^^. 

For any future directed null geodesic k'^ define the quantity ip = —ak°'Ua- Then it 
follows that fc" = —{'ip/a)i°- where 

—r-p / a \ °- 

(14) 

where e = ±1 and the sign of e depends on whether k'^ is outgoing or ingoing. Note 
from the definition of ip that 

a = t± (15) 

where ■?/'+ is the corresponding to n"^ and correspondingly for It follows from the 
form of the LTB metric that 

k'^k'VaUf, = ^ (16) 

Then using the geodesic equation we find 

fc"Va^ = -J^(-r'd + af' + eyTT7a') (17) 

from which it follows that 

rValnV' = -^{r[aa' - aa'] + eJl + fa') (18) 
ar' * 



we will use equation (U81) to find the velocity v of the dust relative to the CMB 
as follows: note that since we have chosen to = it follows that at early times the 
spacetimes becomes FRW. Note also that in FRW spacetimes a' = and therefore ip is 
constant. For CMB photons with the average energy this constant, which we will call 
ijjQ is the same for outgoing and ingoing photons. Integrating equation ( fTSl) back from 
the point at which we want to calculate v to the big bang yields hi{ip/ipo) and since 
we do this for both outgoing and ingoing geodesies, we obtain both ip+/ipo and ip-fipo 
which in turn allows us to calculate and thus the velocity v. 

We will numerically integrate equation (ITSl) . However, in order to do that we will 
need to be able to evaluate all the quantities on the right hand side of that equation 
at all points in the integration. As we will show, all that is needed is to know r and a 
and the rest of the quantities follow. However, r and a also change as the integration 
procedes, so we need to know how they change. From equation ( fT^ we find 



\/ar = eJl + f-r (19) 



rV,a = e^l + /--d (20) 

Note that an ingoing geodesic that goes through r = becomes an outgoing geodesic, 
where this change is implemented as a change in the sign of e. We must evolve the 
system of equations (11811201) . To do this, we must be able to evaluate all quantities on 
the right hand sides of these equations from just r and a, that is, we must be able to 
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find r, /, d,r,a', r', a' and r'. Given r and a we obtain r = r/a and thus we find all 
specified functions of f such as A and /. From equation f|T2l) we obtain 



Ia + - (21) 
a 



(23) 



which then yields r since f = ar. Differentiating equation (1131) with respect to r and 
solving for a' yields 

a' = dIA' (22) 

where 

5 r /„^^a\ 3 /1 + 
'=2A^^[' l' + ^j-2^HT^., 

with s = ^/A/d. This then yields r' through r' = a + fa'. Finally, differentiating 
equation (12T1) with respect to f yields 

, A' /I BI\ , , 

which in turn yields f' through r' = d + fd'. 

What remains is then simply to find the initial values of r and a for each point in 
our past light cone, along with the coresponding values of the redshift z. As shown in 
|T2] this can be done by integrating equations fll9ti20p from our current position, along 
with the evolution equation for the redshift 

rV,^=^(l + z) (25) 
For this integration, the initial value of r is zero, while the initial value of a is given by 

a = \{<:r^l-\) (26) 
where is the ratio of the density in matter to the critical density. 



3. results 

In the models of [12] a choice of '^m is made, and the function A takes the form 

A = — ^ (27) 

l + (cf)' ^ ' 

where the constant c is chosen for best fit with the supernova data. In particular, the 
two models considered are (i) Qm = 0.3, which gives rise to c = 8.5, and (ii) Qm = 0.2, 
which gives rise to c = 5.1. Each of these models fits the supernova data about as well 
as the standard concordance cosmology. Because A — > at large f, these models are 
asymptotically f2 = 1 FRW cosmologies with an underdensity in a region near the center. 
Note that the large value of c means that in these models we live in an underdensity of 
fairly small size. 

Using the method of the previous section, we calculate the velocity of the dust with 
respect to the CMB for each point on our past light cone. The results are shown in 
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Figure 1. Plot of velocity (in units of km/s) vs redshift for the f^M = 0.3 LTB model 
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Figure 2. Plot of velocity (in units of km/s) vs redshift for the Om = 0.2 LTB model 



The motion of galaxy clusters in inhoniogeneous cosmologies 



7 



figures (1) and (2). Here figure (1) is for tlie VLm = 0.3 model and figure (2) is for the 
^A/ = 0.2 model. Also plotted on the figures are the results of observations of the kSZ 
effect for galaxy clusters. [iTl [TBI [I9j 

The observations clearly rule out the VLm = 0.2 model. In addition, the VLm = 0.3 
model does not provide a good fit to the observations, so this model may also be ruled 
out. 

4. conclusions 

Any successful LTB model must pass both the supernova test and the kSZ test. Together 
these tests provide a powerful constraint that rules out the models of [T2]. However, 
the models of [I2] were chosen on the basis of simplicity and fitting the supernova test; 
so it is not surprising that one additional test suffices to rule them out. It is therefore 
possible that reasonable LTB models could be made that have fair agreement with 
both the supernova data and the kSZ data. However, as kSZ measurements improve in 
accuracy, this test will become more stringent and might eventually rule out all LTB 
models. 
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